Turnstile pumping through an open quantum wire 
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We use a non-Markovian generalized master equation (GME) to describe the time-dependent 
charge transfer through a parabohcally confined quantum wire of a finite length coupled to semi- 
infinite quasi two-dimensional leads. The quantum wire and the leads are in a perpendicular external 
magnetic field. The contacts to the left and right leads depend on time and are kept out of phase 
to model a quantum turnstile of finite size. The effects of the driving period of the turnstile, 
the external magnetic field, the character of the contacts, and the chemical potential bias on the 
effectiveness of the charge transfer of the turnstile are examined, both in the absence and in the 
presence of the magnetic field. The interplay between the strength of the coupling and the strength 
of the magnetic field is also discussed. We observe how the edge states created in the presence of 
the magnetic field contribute to the pumped charge. 

PACS numbers: 73.63.Nm 73.23.Hk 85.35.Be 



I. INTRODUCTION 

The time-dependent properties of semiconductor 
nanostructures and their response to electric pulses are 
currently being studied through transient current mea- 
surements and in a pump-and-probe configuration.^^ 
Along with these experimental developments theoreti- 
cal schemes for the description of time-dependent trans- 
port emerged. The methods include the non-equilibrium 
Keldysh-Green function formalismj^ scattering theory,'^ 
and more recently the generalized master equation 
(GME) adapted for electronic transport .I^K^ 

These methods were used mostly for studying the tran- 
sient currents generated by a time-dependent potential 
applied on the sample or by a time-dependent coupling 
between the leads and the sample. An example in the 
first category is the time dependent pumping of electrons 
through a small open system.^ In the second category we 
mention the transient currents and the geometrical effects 
imposed by the lateral confinement of the sample .'^^ In 
another recent work we investigated the modulation of 
the drain current when a sequence of square pulses is ap- 
plied to the source probe connected to a quantum dot and 
a short quantum wire described within lattice model.^ 
That study was motivated by the experiments of Naser 
et aP 

In the present work we further exploit the GME 
method and study the transport properties of a quan- 
tum wire operating in a turnstile regime. The turnstile 
pump is a single-electron device where the sample is pe- 
riodically connected and disconnected with the left and 
right lead respectively, but with a relative phase shift. It 
was experimentally created by Kouwenhoven et al}'^ by 
modulating in time the two tunneling barriers between a 
quantum dot and two leads. The electrons were driven 
by a finite bias between the leads. This setup is different 
from a quantum pump where a current is generated by 
asymmetric external oscillations, but without a bias. In 



the experiment of Kouwenhoven et al. the barrier heights 
oscillate out of phase in the following sense: On the first 
half-cycle electrons enter from the source probe in the 
system but there is no current in the drain probe be- 
cause the corresponding tunneling barrier is high enough 
to prevent this. In the second half-cycle the source is 
disconnected, the drain contact opens, and a discharge 
of the dot follows. It was found that an integer number 
of electrons are transmitted through the structure in each 
pumping cycle. 

More recently, due to the general interest in applica- 
tions of nanoelectronic devices, more complex turnstile 
pumps have been studied by numerical simulations, like 
one-dimensional arrays of junctionals or two-dimensional 
multidot systems. In the present paper we predict that 
the turnstile operation can also be performed in a quan- 
tum wire sample and in an external magnetic field as 
well. Our results are obtained using a parabolic lat- 
eral confinement model both for the quantum wire and 
for the leads. We put a special effort on describing 
the lead-sample contacts. In a previous work we stud- 
ied the turnstile transport through a sample described 
by a lattice (tight-binding) model. The sample was a 
one-dimensional system with two or three sites and the 
transport calculations were done using nonequilibrium 
Keldysh-Green functions. The quantum wire consid- 
ered in this work is much more complex. We start from 
the single-particle Hamiltonian of a two-dimensional wire 
of length Lx parabohcally confined along the y direction 
and with hard- wall conditions at The eigenfunc- 

tions of the Hamiltonian were described in detail in Ref. 
6 and will not be repeated here. 

The material is organized as follows: In Section [ll| we 
briefiy review the main equations of the model and of 



the GME method. Section III is devoted to the numerical 



results for zero magnetic field (III A), in the presence of 



a magnetic field (IIIB), and to the edge states (IIIC). 



The conclusion are given in Section IV 
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II. THE MODEL 

We consider an isolated finite quantum wire of length 
Lx = 300 nm, extended in the x-direction. The width 
of the wire is defined by a parabolic confinement po- 
tential in the ^-direction with the characteristic energy 
HVLq = 1.0 meV. The quantum wire is terminated at 
±Laj/2 with hard wall potentials. In addition, we have 
two semi-infinite leads, one extended from —Lx/2 to — oo, 
and the other one from +1/^^/2 to +oo. Both have a 
parabolic confinement in the y direction with an energy 
of 0.8 meV, and are also terminated at ^L^ll with hard 
walls. The leads and the finite quantum wire, or the 
system, are all subjected to an external constant mag- 
netic field B = 5z. The length of the wire, L^^ and the 
magnetic length modified by the parabolic confinement 
a^y = ^Jh| {m*Qw), with = and the cyclotron 

frequency Uc = eB/{m*c), are convenient length scales 
in the calculations. We assume the GaAs effective mass, 
m* = O.OGTmg. The many-electron Hamiltonian of the 
system composed by the semi-infinite leads and the finite 
quantum wire, but isolated from each other, is: 



(1) 



g,Z=L,R 



where an electron in the system is created (annihilated) 
by the operators (d), and in the leads by (c). Ea 
are the energies of the single-electron states labeled with 
a = 1, 2, 3, ... in increasing order, ei{q) is the energy spec- 
trum of the left and right leads, labeled as / = L and 
/ = R respectively, and q representing a discrete label 
of subbands and a continuous quantum number labeling 
states within each subband. At t = the system is cou- 
pled to the leads with the Hamiltonian 



I 



q,a 
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{TUdic^i}, (2) 



with Xi{t) describing the time-dependence of the cou- 
pling, such that xi{t < 0) = 0^ describing the 
coupling strength of state a and q in the system and the 
leads, respectively. 

The coupling is defined by a nonlocal overlap integral 
of the two states in the region of contact around 
The coupling coefficients are defined phenomenologically 
with the tensoifl^ 



the tunneling rate at the contact between each lead and 
the sample, and also by the parameters 82 and 
which adjust the spatial overlap of lead and sample wave 
functions in the contact region. Since all we have from 
our model are the wave functions derived for the uncou- 
pled subsystems, this intuitive ansatz, Eqs. (|3| and Q, 
is a convenient way of describing the coupling. In our 
calculations these coupling parameters will be the same 
for both leads and hence the label / will be omitted. 

The system can be subjected to a bias A/i = ijll — I^r 
and in order to reduce the number of many-electron 
states (MES) to a reasonable number in the follow- 
ing calculations we limit the number of single-electron 
states (SES) for the particular calculation by selecting a 
window of relevant states around the bias window, i.e. 
[11 R — A, /ii, + A], such that the transport properties are 
not changed significantly by extending the window. We 
consider these states relevant for the transport, or "ac- 
tive" , while all the other states are "frozen" , being either 
permanently occupied or permanently empty. In addi- 
tion to the electrons frozen in the states below the active 
window (which are not included in the transport calcu- 
lation), we also assume the lowest active state occupied 
at the moment t = i. e. when the contacts begin to 
operate. In this way the transient phase is shorter than 
if we would assume an empty active window, and we can 
spend less computing time until the system reaches the 
periodic state. 

The time-evolution of the total system - finite wire 
and leads - after the coupling at t = 0, can be described 
by the Liouville-von Neumann equation for the statis- 
tical operator W{t). The evolution of the finite wire 
itself can be captured by the reduced density operator 
p{t) = Tii^TiYiW (t) (RDO), i. e. by averaging over the 
the lead variables. From the resulting integro-differential 
equation we retain only the lowest order (quadratic) 
terms in Ht and obtain^^ 



- ^ E fdqXim[Ti,n,iit)] + h.c.), (5) 

where we have introduced two operators to compactify 
the notation 



TL,= [ rfr.r'(*^(r'))**f(r),i,(r,r') + /..c., (3) ^^t) = UlU) [^ds Ms)U,^is) 



Ms-t)/h)Ei{q) 



where a nonlocal-overlap of the wave functions in the 
system and the leads is modeled by 

5i,(r,r') = 5^exp [-5[{x - x'f - 5[{y - y'f] 
\Ea-e\q)\\ 



n,;(s) = Us{s) (7;V(s)(l - fi) - Pis)T^fi) ulit), 

with Us{t) = e^'^^/^^^s^ Q^Yid a scattering operator T act- 
ing in the many-electron Fock space of the system 



exp 



(4) 



a, (3 



The strength of the coupling between the leads and the 
sample is defined by the parameter g^^ which captures 



(6) 
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With the RDO it is possible to calculate the statistical 
average of the charge operator Qs = e'^^d^da for the 
coupled system 



{Qsit)) 



Tv{W{t)Qs} 
Tvs{p{t)Qs} ^ 



Trs{[TrLRl^(t)]Qs} 
eJ2i^A,,\p{m, (7) 



with the traces assumed over the Fock space. The average 
time-dependent spatial distribution of the charge can also 
be obtained, 

(Qs(r,i)) = e^^*:(r)*6(r)p^,(t)(z.|44|M)- (8) 

ab fJ^i^ 



The net current flowing into the sample is 



{j{t)) 



(Mt)) - { JrW) 

d{Qs{t)) 



(9) 



dt 



The total current in Eq. (9| is given by the left hand side 
of the GME, Eq. ([5|, and the partial currents associated 
to each SES and each lead correspond to the terms of the 
sum on the right hand side (the trace of the commutator 
of p and Hs is zero). 

The functions which modulate the coupling between 
the leads and the quantum wire are built in the following 
way: For < t < Ti we use the Fermi-like function 
f{t) = (e^^ + with 7 = 1.0 ps~-^, and we define 

Xi{t) = 1 — where / = L, R is the lead index. Then, 

for t > Ti^ Xi{t) become step-like functions alternating 
between and 1, both with the same period T, but with 
a delay of T/2. In this way we mimic the on/off contact 
switching done in the turnstile experiments. We choose 
< = + ^/2, which means we first switch off the 
left contact while the right one is still on. Then the left 
is turned on again while the right is turned off, and so 



III. NUMERICAL CALCULATIONS, RESULTS, 
AND DISCUSSION 



A. No magnetic field 

The energy spectrum of the leads and of the sample are 
shown in Fig.[T] The leads have parabolic subbands while 
the sample has discrete levels. The maximum energy for 
each subband shown in the graph indicates the corre- 
sponding maximum wave vector in the (^a^^-integration of 
the GME. The chemical potentials in the leads defining 
the bias window (BW) are shown with the dashed hori- 
zontal lines. We consider two BW's: BWl with chemical 
potentials jul = 1-48 meV and pR = 0.78 meV, and BW2 
with fiL = 2.48 meV and fiR = 1.78 meV respectively. 
In both cases the applied bias is eVbias = Pl — I^r = ^-"^^ 
meV. In the numerical calculations of the reduced sta- 
tistical operator we also include the sample states with 



energy outside the BW, between the limits /j^r — A and 
/iL + A with A = 0.1 meV. The first active window con- 
tains 4 SESs and the second one contains 5 SESs. 
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FIG. 1: (Color online) The energy spectrum of the leads 
vs. the scaled wave vector qaw The subband index is 
n = 0,1,2,3. With dots the energy spectrum of the iso- 
lated sample. With horizontal wide dashed lines the chemical 
potentials ijll — 1.48 meV and iir — 0.78 (BWl), and with 
horizontal narrow dashed lines fiL = 2.48 meV and fiR = 1.78 
meV (BW2). 

In Fig.|2]we show the time-dependent total occupation 
nt of the relevant (active) SES for BWl and BW2. The 
electrons occupying the states situated below the active 
windows are not included, but only those on the relevant 
states (or within the active windows). In this example 
the parameters characterizing the coupling of the sample 
to the leads are g^aw = 926 meV (with a^j measured in 

nm), Sio^ = 1.0, and 620^ = 2.0. The large value for 

3/2 

goaj is consistent with the small overlap of the wave 
functions from the lead and from the sample mandated 
by the large values for ^la^ and ^2<^^- The largest con- 
tribution to the overlap comes from the "contact" area 
within one around the ends of the finite quantum 
wire at X = iL^^. The unusual dimension of go^ which is 
energy x length"^^^ is a result of different dimensions 
of the wave functions in the lead and in the sample: the 
former is length ^ , being unbounded in the direction 
along the lead, and the later is length"^. We also select 
A^; = 0.25hClw in Eq. Q. The time-dependent switching 
functions Xi{'^) also shown in the figure, with period 
T = 60 ps, for / = L, R. 

Comparing the charge oscillation for the first and for 
the second active window we see that more charge is 
transferred through the system for BW2 than for BWl. 
There are two reasons for that. One reason is that BW2 
includes three subbands of the leads, whereas BWl in- 
cludes only two (Fig. [T]). The second reason is that the 
energy dispersion in each subband increases with increas- 
ing energy, such that the electrons have higher speed in 
BW2 than in BWl, and thus an increased contribution 
to conduction. 

It is interesting to observe the behavior of the states 
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FIG. 2: (Color online) The time dependent total occupation 
of the relevant SES for a system occupied initially with 1 
electron (/io = 2) for the first (solid), and second (dotted) bias 
window. The time coupling functions Xi{^)i with I — L^R^ are 



shown for reference. Other parameters: g^d 



3/2 



926 meV, 



^iCiw = 1-0, fca^ = 2.0, and T = 60 ps. In each case there is 
one electron on the lowest relevant state at the initial moment. 



situated at the boundaries of the BW. To show that we 
choose the BWl and display in Fig. [3^ the partial cur- 
rents created by each sample state included in the calcu- 
lation. These are the states 2,3,4,5 in the order of the en- 
ergy, shown with dots in Fig.jl] State 1 (the lowest dot) is 
considered totally occupied and frozen. In this case state 
2 is slightly below the window, state 3 is within the win- 
dow, and states 4 and 5 are slightly above the window. 
In Fig. we show the same partial currents, but now 
for a slightly higher chemical potential of the left lead, 
jiL = 1-54 meV, instead of /i^ = 1-48 meV used in Fig. 
[3^. The state number 4 is now included in the BW and 
consequently the corresponding current increases. The 
current of the state 5 also increases a little bit, while the 
current associated to the other states does not change. A 
similar behavior is displayed by state number 9 situated 
on top of BW2 (not shown). 

The efficiency of the turnstile operation depends on 
the pulse length T. The previous results are obtained 
with T = 60 ps. In those setups the system transfers at 
least two electrons per cycle. The present GME method 
is valid in the lowest (quadratic) order of the lead-sample 
coupling, which means the tunneling of the electrons from 
the leads to the sample and back is a relatively slow pro- 
cess. Therefore by increasing or decreasing the pulse du- 
ration the transferred charge increases or decreases re- 
spectively. Denoting by Tt the characteristic tunneling 
time, if T < the turnstile operation is not expected to 
work, the allowed time for the charging and discharging 
of the system being too short. This is the case for T = 10 
ps as shown in Fig. [4j when clearly very little charge can 
enter and leave the system in a pumping cycle. Actually, 
the time Tf depends on the pairwise coupling (or overlap) 
of each state of the leads to each state of the sample with 
energies within the BW. But in order to obtain signifi- 
cant pumping effects, the pulse duration has to include 
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FIG. 3: (Color online) The partial current entering the sys- 
tem from the left lead when the bias window contains one 



respectively two SES. Other parameters: god 
Siai = 1.0, 62ai = 2.0, and T = 60 ps. 
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926 meV, 



the time of flight (or propagation time) of electrons along 
the wire. This extra time depends on the energy of the 
electrons injected from the left lead all the way into the 
right lead. 

Therefore, for a longer pulse, T = 40 ps, the turnstile 
pumping process is able to transfer charge through the 
sample. The occupation number has a triangular shape 
in time and it becomes periodic after 2-3 cycles. Dur- 
ing the initial cycle, which includes the initial charging 
phase, the system accumulates more than 2 electrons and 
the steady state is already reached at about 30 ps when 
the charge in the system is saturated. This is possible be- 
cause the right contact is still off. The right contact opens 
for the first time at 45 ps allowing more than 1 electron 
charge to pass into the right lead. For a longer period, 
like T = 60 ps, the occupation number develops toward a 
saw-tooth profile typical for the charging/relaxation pro- 
cesses. The asymmetry of the charge peaks is determined 
by the direction of the bias: the electrons leave the sam- 
ple faster than they entered. The system drives almost 
two electrons from one lead into the other, which is re- 
markable given the length of our sample (300 nm). 

Another important aspect in our model is the strength 
of the lead-sample coupling, i. e. the parameters go^ Ji, 
and 62 in Eq. (pi). Our GME implementation is restricted 
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FIG. 4: (Color online) The time dependent total occupation of 
the sample for three pulses at zero magnetic field. The trace 
for T = 60 ps is the same as in Fig. [2] Other parameters: 



3/2 



: 926 meV, 5iai = 1.0, and 52ai = 2.0. 



ence of a magnetic field shown in the next examples, the 
scaled parameters goctw'^ and ^i,2<^^ chosen such that 
the physical values go and ^1^2 are the same as for B = 0.9 
T. The time dependent occupation of the states in the 
active window is shown for longer times than in the pre- 
vious figures to indicate better the final periodic regime. 
It is not surprising to see that the pumping amplitude 
increases with the coupling strength, since tunneling of 
electrons becomes more likely. The same can be said 
about the current, which is essentially the time deriva- 
tive of the charge, and it is visible in Fig. [5] that the 
slope of the charge increases with the coupling strength. 
The increase of the transient current with the coupling 
strength has also been shown recently by Sasaoka et al. 
although in a quantum dot pumping device.^ 



B. Magnetic field present 
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FIG. 5: (Color online) The time dependent total occupation 
for 5 = T and two values of the coupling strength. Weak: 
goot/^ = 1408 meV, ^la^ = 1.85, ^2^^ = 3.7. Intermediate: 
goal^^ = 1408 meV, 6ial = 1.39, 62ai = 2.77. Strong: 
goaK'^ = 1824 meV, 6iai = 1.39, 62ai = 2.77. The pulse 
period is T = 60 ps. 



to the lowest order in Ht (quadratic), supplying the 
integro-differential equation ([5|, and thus the parameters 
have to be appropriately selected. In general it is difficult 
to evaluate whether the coupling strength is sufficiently 
low. A necessary (although not sufficient) condition is to 
obtain positive diagonal elements of the statistical oper- 
ator, which are the populations of the MES and hence 
probabilities. Although we always check in our calcula- 
tions, strictly speaking this condition does not guarantee 
the validity of the lowest order approximation (quadratic 
in Ht)- So in practice we cannot avoid choosing our pa- 
rameters in a semi-empirical manner. To have an idea 
about the relation between the pumping amplitude and 
the coupling strength we show in Fig.|5]three calculations 
done with three strengths of the coupling, which we con- 
sider in relative terms weak, intermediate, and strong. In 
order to compare with the results obtained in the pres- 



Following our ansatz used to describe the system-leads 
coupling, Eqs. ([3|-(|4|, we can visualize systems where the 
physical parameters go^ (^1, and S2 should be assumed 
constant, and others where the scaled values goa^'^ ^ 
(5ia^, and (^2<^^ should be kept constant with changing 
values of the magnetic field rather. In the following we 
explore both possibilities. Our choice of A^; = 0.2bh^}w 
depends also on the value of the magnetic field, but we 
have checked that variations of this parameter between 
a scaled version and a fixed physical one will only lead 
to vanishing quantitative changes in calculations where 
we consider a range in the q integration of the GME ([5| 
that includes 4 subbands of the leads. This choice of A^; 
guarantees always the same number of subbands in the 
calculation. 

In the presence of a magnetic field, in a first approx- 
imation one might expect the amplitude of the charge 
oscillations to decrease. The reason being that the elec- 
tronic trajectories bend due to the Lorentz force and the 
electrons might return to the source lead rather than 
traveling directly to the drain lead. At the same time we 
know that a magnetic field generally reduces backscatter- 
ing. In Fig. [6] where we show the time dependent total 
charge for different values of the magnetic field we see a 
reduction of the charge oscillations with increasing mag- 
netic field. For comparison we include the B = case, 
also shown in Fig. [2] and Fig. [4] To compare the re- 
sults with and without magnetic field we initially keep 

the same coupling parameters as measured with respect 

3/2 

to the modified magnetic length, a^^ i. e. godw = 926, 
meV Sia^^ = 1.0, and S2a'^ = 2.0. The modified mag- 
netic length decreases with increasing magnetic field, 
as do the features of the wave functions. We select the 
bias windows like BW2, with three states included in the 
window and two marginal states in the extended regions. 
The energy levels and the wave functions depend on the 
strength of the magnetic field both in the sample and in 
the leads. In order to make the comparison shown in Fig. 
|6] meaningful for each value of the magnetic field we shift 
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the chemical potentials in the leads such that the same 
energy levels of the quantum wire sample are contained 
in the bias window with a fixed width. 
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FIG. 6: (Color online) The time dependent total occupation 
of the relevant SESs for different values of the magnetic field. 
go(hJ^ — 926 meV, SiO^ — 1.0, 620^, — 2.0, and T = 60 ps. 



4.2 
4 
3.8 




40 80 120 160 200 240 280 320 360 400 
t(ps) 

FIG. 7: (Color online) The time dependent total occupation 
for B = 0.9T and three coupling strength: weak {goaK'^ = 926 
meV, 6iai, = 1.0, and 620^, = 2.0, intermediate {qqo^'^ — 926 
meV, 5ia^ — 0.75, and ^2^^ — 1-5), and strong (goaw'^ — 
1200, 5ia^ — 0.75, and fca^ = 1.5). The pulse period is 
T = 60 ps. 



The efficiency of the turnstile pump tends to decrease 
with increasing the magnetic field. For 5 = an amount 
of charge Qp ^ 1.8 electron can be transmitted along 
the wire sample in one cycle, and for 5 = 0.2 T 
1.2. For stronger magnetic fields Qp drops to 0.8 for 
B = 0.4 T, and Qp < 0.5 at 5 > 0.6 T. It is also clear 
that, at least for the present coupling parameters, the 
charging time increases in the presence of the magnetic 
field. During the initial cycle, i. e. for t < 60 ps, there 
is less charge accumulated in the system in the normal 
switching regime than for 5 = 0, and also the charging 
process continues even after the pumping begins. The 
periodic regime is reached after a number of cycles which 
increases with B. 

Although the presence of the magnetic field reduces the 
electron transfer, the pumped charge still increases with 
increasing coupling between the leads and the sample. 
To show that we solve the GME for a fixed magnetic 
field B = 0.9 T for three sets of coupling parameters, 
which we again call (in relative terms) weak, intermedi- 
ate, and strong coupling, respectively, see Fig. [T) The 
parameters for the weak coupling are the same as in Fig. 
[6j For the intermediate coupling we use g^aw^ = 926 
meV, ^la^ = 0.75, and S2a'^ = 1.5. For the strong cou- 
pling goctw'^ = 1200 meV, 5ia2 = 0.75, and (^2<^^ — l-^- 
These parameters correspond to the same physical pa- 
rameters that were used in Fig. [5] So, if we now compare 
the results at 5 = 0.9 T with the results at 5 = we 
see quite similar charge amplitudes, but somewhat more 
sensitive to the coupling strength at B = 0.9 T. For ex- 
ample, at 5 = 0.9 T we obtain Qp ^ 0.4 electrons at 
low coupling, Qp ^ 1.2 at intermediate coupling, and 
Qp ^ 2.4 at strong coupling. For B = these numbers 
are Qp ~ 0.3, 1.2, and 1.8, as can be read from Fig. [5] 

We believe the increased sensitivity of the pumping to 
the contact strength in the presence of the magnetic field 



is a result of the edge states created in the sample. The 
edge states are indeed more and more pronounced with 
increasing magnetic field, and so is the pumped charge 
if the edge states are in good overlap with the coupling 
region. Therefore with magnetic fields of about 1 T and 
strong coupling we can obtain Qp of about 2 electrons. 
Fig. [8| which is slightly more than at zero magnetic field 
but with weak coupling Qp ^ 1.75, visible in Fig. [7[ 

The time dependent total currents in the left and right 
leads are shown in the lower panel of Fig. [S] for 5 = 1 
T. The currents suddenly vanish at each contact when 
the contact is closed. The currents have a trapezoidal 
shape for the pulse duration t = 60 ps, but they become 
triangular for shorter pulses like T = 10 ps (not shown). 



C. Analysis of edge states 

Clearly, states in the sample with higher probability 
at the contact edges contribute more to the pumping. In 
Fig. |9] we show the coupling strength of each of the five 
states of the finite wire involved in the transport for each 
magnetic field of Fig. [8) only for the lowest subband of the 
leads. The general trend of the coupling coefficients is to 
decrease with increasing magnetic ffeld. So the decrease 
of the pumping when the magnetic ffeld increases can be 
explained by the decrease of the coupling strength. The 
contribution of each state to the transport is given by an 
integration over all the lead states q in the GME. 

To see the contribution of each state to the transport 
we show in Fig. [lO] the time dependent occupation of all 
5 states included in the calculation for 5 = 1.2 T, which 
are the states number 6-10 in the single-electron energy 
spectrum of the ffnite wire. It is clear that the three mid- 
dle states a = 7, 8, 9 contribute unequally to the pump- 
ing. These states are well inside the BW, but the coupling 
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FIG. 8: (Color online) (a) The time dependent total occupa- 
tion of the relevant SESs for B = 0.9-1.2 T. (b) The total cur- 
rent entering the system from the left lead and the total cur- 
rent exiting the system into the right lead for B = 1.0 T. Pa- 



(b) 



Left 
Right 



rameters: goalC^ = 1200 meV, 6iai 
The pulse period T = 60 ps. 



0.75, and S2a^ = 1.5. 



energies are slightly different as seen in Fig. |9] for B = 1.2 
T. For the other values of the magnetic field shown in Fig. 
[9] the states within the BW have nearly equal coupling 
to the leads and consequently nearly equal contributions 
to transport. So in general the coupling strength may 
depend on the state and so does the corresponding par- 
tial current. The states 6 and 10 included in Fig. [lO] are 



slightly outside the BW and obviously their contribution 
to transport is smaller. 

In Fig. [11] we show the probability density associated 
to the five active single-electron states (number 6-10) of 
the finite quantum wire. The figures indicates that all 
the five states have the characteristics of an edge state. 
In vanishing magnetic field the probability density of the 
active states is on the average smeared over the finite 
wire. As the magnetic field increases the Lorentz force 
squeezes the probability of some states close to the edge 
of the finite wire. This also happens at the hard- wall ends 
of the wire, the contact area. This fact explains why the 
increasing of the coupling through increasing should 
be more effective at high magnetic field. From Fig. pTjit 
is evident that the edge states will have different coupling 
strengths to the leads due to the difference in their finer 





FIG. 9: (Color online) The coupling energies between the 
states a = 6 — 10 and states of the lowest subband in the left 



(or the right) lead for B = 0.9 - 1.2 T. ^0- 
Sia'i, — 0.75, and ^2^^ — 1.5. 
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FIG. 10: (Color online) The time dependent total occupation 



of the relevant SESs for B = 1.2 T. goa 
^la^ — 0.75, and ^2^^ — 1.5. T = 60 ps 
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1200 meV, 



structure in the contact area. This finer structure in the 
contact area of the wire induces differential coupling to 
the states in the different subbands of the leads. 

The time-dependent charge in the quantum wire is 
shown in Fig. [12] for the parameters used in Fig. ^jp. 
The charge distribution reflects the geometry of the quan- 
tum wire and of the five SESs involved, and indicates the 
propagation of the electrons in the system. The selected 
time moments cover the initial charging cycle plus a part 
of the next cycle. It is interesting to observe how the 
electrons are injected at the left contact into the sam- 
ple traveling along the quantum wire on the upper edge 
channel, and how they are reflected back at the right 
contact traveling along the lower channel. 



8 





FIG. 11: (Color online) The probability density of the single- " " " X„ ' ' " " " 1. ' ' " " " x/L ' ' 

electron eigenstates of the sample number 6-10 (top left - right 

- down), for B = 1 T. FIG. 12: (Color onhne) The average spatial charge distribu- 

tion in the quantum wire for B = 1 T with strong coupling, 
at different moments of time. 

IV. CONCLUSIONS 



The turnstile conduction of a two-dimensional 
parabolic wire, seen as "the sample", attached to semi- 
infinite leads of a similar parabolic shape has been stud- 
ied using a non-Markovian generalized master equation 
method. This system is far more complex than the 
simple two- or three-level system with one-dimensional 
leads considered in an earlier publication based on the 
Keldysh-Green functions approach. The eigenstates of 
the leads and of the finite wire have been calculated in 
an external perpendicular magnetic field using a com- 
bination of analytical and numerical methods for large 
functional basis sets. We have taken into account the 
subband structure of the leads to which it is attached. 
We have also described phenomenologically the coupling 
coupling between the states in the leads and the states in 
the finite wire as a nonlocal overlap of the wave functions 
from both sides of the contact. 

We have analyzed the effects of the bias window, pulse 
length, and magnetic field on the evolution in time of the 
number of electrons in the sample. We have found that 
longer pulses are more favorable for the turnstile pump- 
ing as the electrons need time to propagate along the 
sample wire. One or two electrons could be transferred 
through the 300 nm wire using pulses of 40 or 60 ps. 

The comparison of the results obtained with and with- 
out a magnetic field is a difficult issue. If a magnetic 
field is present all energies shift, all wave functions change 
(both in the sample and in the leads), and also all ele- 
ments of the coupling tensor T^^ between lead and sample 
states change. Of course these changes depend on the 



strength of the field. Therefore it is difficult to create 
similar conditions for two different field values, i. e. the 
same number of states in the bias window and the same 
coupling energies, in order to compare only the ampli- 
tude of the pumped charge. To do that we have selected 
the parameters describing the phenomenological coupling 
of the leads to the finite wire in two different ways, both 
scaled and not scaled with the effective width of the sam- 
ple, which implicitly depends on the magnetic field. This 
is an issue that can only be better resolved with a more 
involved microscopic model of the coupling and compari- 
son to experiments where the coupling strength could be 
varied, for example by using finger-strip gates. 

The charge distribution inside the system. Fig. 
[T2j emphasizes the dynamics induced by the charg- 
ing/discharging sequences. The charge propagation 
along edge-states in a strong magnetic field indicates 
that the optimal turnstile period depends on the external 
magnetic field. Experimental studies of turnstile pump- 
ing in quantum wires have to clarify the relationship be- 
tween the magnetic field, pumping amplitude, and con- 
tact strength. 
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